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Abstract 


Finite element method (FEM) has been a very powerful tool to solve 
many complex problems in electromagnetics. The goal of the current 
research at the Langley Research Center is to develop a combined 
FEM /method of moments approach to three-dimensional scattering/ 
radiation problem for objects with arbitrary shape and filled with complex 
materials. As a first step toward that goal , an exercise is taken to 
establish the power of FEM, through closed boundary problems. This 
paper demonstrates the development of FEM tools for two - and three- 
dimensional eigenvalue problems in electromagnetics. In section 2 , 
both the scalar and vector finite elements have been used for various 
waveguide problems to demonstrate the flexibility of FEM. In section 3, 
vector finite element method has been extended to three-dimensional 
eigenvalue problems . 


1. Introduction 

The finite element method (FEM) has been 
widely used as an analysis and design tool in many 
engineering disciplines like structures and computa- 
tional fluid mechanics. Though FEM has been ap- 
plied to electromagnetic problems, it was mainly con- 
fined to electrical machines and magnetics (ref. 1). 
In the past 20 years there has been a great inter- 
est in application of this method to microwave com- 
ponents such as waveguides and antennas. But for 
many years, its use has been restricted because of the 
so-called spurious solutions in vector finite elements 
(ref. 2). Very recently, the “edge elements” have been 
employed successfully for vector formulations with- 
out resulting in “spurious solutions.” In the recent 
past, use of these edge elements in FEM has revived 
an interest in applying FEM to microwave engineer- 
ing problems (ref. 3). This in combination with the 
advances in computer hardware and software helped 
to make FEM an attractive tool for electromagnet- 
ics. Also, there are a variety of commercial geomet- 
rical modelling tools to accurately model any three- 
dimensional geometry and to generate the required 
mesh with any kind of elements such as triangles and 
tetrahedrals (refs. 4 and 5). 

In this paper, the FEM tools for analyzing eigen- 
value problems in electromagnetics have been de- 
scribed. This paper is divided into two parts: sec- 
tion 2 deals with the two-dimensional problems; 
section 3, with the three-dimensional problems. 
Throughout this paper triangular elements are used 
for modelling two-dimensional problems and tetra- 
hedrals are used to model the three-dimensional 
problems. 

In section 2.1, a scalar FEM formulation is used 
for two-dimensional arbitrarily shaped waveguides. 
Triangular elements with nodal basis functions are 


used to formulate the FEM matrices. The eigen- 
values for different types of waveguides are obtained 
and the field intensity plots are presented for various 
waveguide modes. 

In section 2.2, a vector FEM is introduced with 
two-dimensional edge elements for analyzing in- 
homogeneous waveguides. For the sake of clarity in 
formulation, section 2.2 is divided into four sections. 
Section 2.2.1 gives the solution of homogeneous wave- 
guide problem with two-component transverse vector 
fields. Section 2.2.2 gives the calculation of eigen- 
values for inhomogeneous waveguides using the three- 
component vector fields. Combination of edge and 
nodal basis functions have been used for transverse 
and longitudinal field components, respectively. Sec- 
tions 2.2.3 and 2.2.4 extend the formulation in sec- 
tion 2.2.2 to determine either the wave number or 
the propagation constant for inhomogeneously filled 
waveguides when one of them is specified. 

In section 3.1, formulation for three-dimensional 
vector FEM is described. Edge basis functions for 
tetrahedral elements are introduced to formulate fi- 
nite element matrices for three-dimensional cavities 
filled with inhomogeneous material. 

In sections 2.1, 2.2, and 3.1, numerical examples 
are presented to show the validity of the analysis and 
the computer programs developed. For all the ex- 
amples, FEM results in good accuracy. All the nu- 
merical examples have been checked for numerical 
convergence. By virtue of FEM, the computer codes 
presented in this paper can handle any arbitrarily 
shaped geometries filled with inhomogeneous mate- 
rials, unless otherwise mentioned. 



2. Two-Dimensional Problems 

2*1. Homogeneous Waveguides — Scalar 
Formulation 

This section deals with the solution of two- 
dimensional waveguide problems with closed bound- 
aries using the Galerkin finite element method. The 
wave equation is solved for a generalized prob- 
lem with nodal-based first-order triangular elements 
(ref. 1). Finite element matrices are derived and a 
computer program to calculate the eigenvalues and 
electric field distributions is presented. 

2.1.1. Formulation 

The scalar potential function ip satisfies the 
Helmholtz equation with wave number k c 

+ kfy = 0 (1) 

within and on the boundary dT indicated in figure 1 . 
This is the “strong” form of the scalar Helmholtz 
equation. In the strong form, the unknown appears 
within a second-order differential operator. To make 
the equation suitable for a numerical solution, it can 
be converted into the “weak” form by multiplying 
both sides with a testing function T s and integrating 
over the surface T; that is, 

J I [t s (vfy) + k 2 T s i>] ds = 0 (2) 

r 

The first term in equation (2) can be written as 

I Jt s (vfy) ds = f f T S (V< • Vti>) ds (3) 
r r 

The following vector identities can be used to modify 
equation (3): 

V* * [T,(V t tfO] = V t T s • + T s (V t • V s V0 (4) 

and 

Vi*V^ds= [vtipridl (5) 

r dr 

where n is the unit normal along the boundary dT. 
Equation (2) can now be written as 

JJ(V t T s -V t ip)ds-k 2 J j T s ij} ds = J T s ^ dl 

T TdT 

( 6 ) 


where d'lp/dn is the normal derivative of \p along 
the boundary dT. The term on the right-hand side 
vanishes as T s vanishes on the PEC boundary for the 
TM case and d'lp/dn vanishes at the boundary for the 
TE case. Hence equation (6) can be written as 

J f (V t T s • V#) ds = k 2 c J J T s il>ds (7) 

r r 

2.1.2. Discretization 

The problem region is discretized with the first- 
order triangular elements. Within the triangular ele- 
ment given in figure 2, 'ip is adequately approximated 
by the expression (ref. 1) 

'ip = a + bx + cy (8) 

The solution is piecewise planar but continuous 
everywhere. At the vertices 1, 2, and 3, the potential 
can be expressed as 

ip\=a + bxi+ cy\ (9) 

xp 2 = a + bx 2 + o/2 (10) 

^3 = a + &X3 + q/3 (11) 



Figure 1. Geometry of problem. 



Figure 2. Single triangular element. 
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IT 1 


From equations (9), (10), and (11), the coefficients 
a, b , and c are evaluated as 


( 12 ) 


Hence, equation (8) can be rewritten by substituting 
for a, 6, and c 


and the right-hand side as 


a 


"1 XI yi' 

-1 

>l' 

b 

= 

1 x 2 1/2 


i >2 

c_ 


.1 x 3 y 3 _ 


> 3 . 



'a' 


"1 X\ yi~ 

-1 

■tfr 

0 =[i * y] 

6 

= [1 1 y] 

1 ^2 1/2 


02 


_c_ 


- 1 ^3 1/3 - 


-■ 03 - 


(13) 

Equation (13) can be written as 
3 

^ = (14) 

i=i 

where ai{x,y) is given by 

<*i{x,y) = ^{di + bix + ay) (i = 1,2,3) (15) 
and a j , , and q are given by 

<*i = x jVk ~ x kVj (16) 

b i = Vj ~ Vk (17) 

Cl — x k — Zj (18) 


where z, j, and k are cyclical; that is ( i = 1, j = 2, 
k = 3), (z = 2, j = 3, & = 1) and (i = 3, j = 1, 
fc = 2), and A is given by 


A = 


1 


XI 

Vl 



^2 

V2 


(19) 

^3 

V3 



(as 

per 

Galerkin’s 

tech- 

u 

= 1 

,2,3) 

(20) 


nique (ref. 1)), 


and the element representation in equation (14), the 
left-hand side of equation (7) can be evaluated over 
a single element as 


II 


{V t T s ■ V t i>)dxdy 


= X>i/ [ (V(Q; • V t<Xj) dx dy (j =1,2,3) (21) 
i=] v 


J J T s ipdxdy = J J atiOtjdxdy (j = 1,2,3) 


A i=l A 

Hence for each element, equation (7) becomes 
3 


( 22 ) 


E ^ J J ( Va * ■ Va ?) dx d y 

i=i v 

= k 2 c ^2 ipi J J otiOtj dx dy (j = 1, 2, 3) (23) 


3 

I 

<=i a 

And this can be written in a matrix form as 

l2 

Peljm = ' 

where 


[SellM = [T e l] [0] 

[S e i] = J J (Va, ■ Vqj) rfxdy 

A 

[ T el] = JJ dxdy 


doti ^ dc*i^ 

Vai = ^ X+ ^ y 

From equation (15), Va* can be written as 


and hence. 


[Vo* * Vet,] = 


Vai = — (6*$ + Ciy) 


Vai ■ Vo i Vqi * V 02 Vqi ■ V 03 
V«2 * Vai Vq 2 • V02 Vo 2 ■ V03 

LVo3*Vai VQ3 - Vq 2 Va 3 -Va3J 


(24) 

(25) 

(26) 

(27) 

(28) 
(29) 


Substituting equation (28) into equation (29) gives 


[^ Va,] = _ 


bj + c\ b\b 2 + C 1 C 2 6163 4- C 1 C 3 
b 2 bi + C 2 C 1 b\- hc$ b 2 b 3 + c 2 C3 

. b 3 bi + C 3 Cj b 3 b2 + C 3 C 2 4- C 3 

(30) 

The matrix [S e i] can be evaluated by using equa- 
tion (30) to obtain 


[Sell ~ 


J J (Vaj * Vaj)dxdy 

r 

j J ctiCXj dx dy 


= A[Vo t • Vqj] ( 31 ) 


[T e l] 


(32) 


3 



The matrix [T e j] has been evaluated by Silvester 
(ref. 6) and is given in a simple form as 



"2 1 
1 2 
1 1 


1 

1 

2 


(33) 


The matrices [S e J and [T e j] are evaluated for each 
element and are assembled over the entire region 
according to the global node numbering to obtain 
a global matrix equation (ref. 1) as follows: 


walls of the waveguide to satisfy the Dirichlet bound- 
ary conditions for the longitudinal electric field com- 
ponent. A very simple way of implementing this is to 
ignore the nodes on PEC while forming finite element 
matrices. This will result in lower order matrices for 
the TM case than those for the TE case. Once the 
scalar potential is obtained, the transverse electric 
fields for the TM modes are given by 

E x = -Z™ (40) 


[SM = * c 2 [TM (34) 


r 7 tm 

Ey ~ ~ Z ° dy 


(41) 


This results in matrices of the order n x n where n 
is total number of nodes. With equations (30), (31), 
and (33), the eigenvalue equation (eq. (34)) is solved 
for by the standard eigenvalue solvers from the 
EISPACK library (refs. 7 and 8). The cutoff wave 

numbers are then given by \Jk%. 

2.1.3. Field Computation From Scalar Potential 

Once the scalar potential is calculated at every 
node, the electric field can be calculated for both the 
TE and TM modes by the following formulation. The 
scalar potential at any point (x, y ) inside a triangular 
element is given by 


3 

^ ( 35 ) 

i= 1 

These scalar potentials can be differentiated with 
respect to x and y to obtain the following expressions: 


dx 


1 

2A 


3 


i= 1 


(36) 


d4> 

dy 


-s?> 


2=1 


(37) 


For the TE modes, the transverse electric field com- 
ponents E x and E y inside an element are given by 


E x = 


dip 

dy 


(38) 


Ey — 


dip 

dx 


(39) 


When obtaining the scalar potential for the TM 
modes, the scalar potential is set to zero on PEC 


where is the characteristic wave impedance for 
the TM mode. 

2.1.4 • Numerical Examples 

A computer code HELM 10 was written to imple- 
ment the formulation presented in section 2.1. The 
flowchart for the implementation of the FEM solu- 
tion is given in figure 3. Numerical examples for the 
rectangular waveguide, the circular waveguide, and 
the coaxial line are given as follows: 

Rectangular waveguide: The cutoff wave numbers 
k c of a rectangular waveguide calculated by using 
HELM10 are presented in table 1 along with ana- 
lytical results from reference 9. The geometry of the 
rectangular waveguide (or/b r = 2) is shown in fig- 
ure 4. The numerical results presented are achieved 
with 400 triangular elements over the cross section 
of the waveguide. The eigenvectors for some of the 
modes have been calculated and the electric fields of 
the corresponding modes are plotted in figure 5. 



Figure 3. Flowchart for FEM solution. 
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Figure 4. Geometry of rectangular waveguide. 


Circular waveguide: The cutoff wave numbers for 
a circular waveguide of unit radius were calculated 
with HELM10 and compared with analytically avail- 


able data from reference 9. (See table 2.) A cross 
section of the circular waveguide is shown in figure 6. 
Two hundred triangular elements have been used to 
model the geometry. The eigenvectors of selected 
modes have been calculated and the electric fields of 
these modes axe plotted in figure 7. 

Coaxial line: The cross section of the coaxial line 
is shown in figure 8. The HELM 10 program is used to 
calculate the cutoff wave numbers and corresponding 
electric field strength of higher order TE and TM 
modes. A triangular mesh with 340 elements has 
been used to model the geometry. Table 3 presents 
the cutoff wave numbers computed for 7 * 2 /n = 4 by 
HELM 10 and the analytically available cutoff wave 
numbers in the literature (ref. 10). For TM modes 
the potential on the inner and outer conductors is set 
to zero. The transverse electric field components are 
calculated and plotted in figure 9. 


Table 1. Cutoff Wave Numbers for Rectangular Waveguide 


TE 

TM 

JzqCL'p 

Analytical (ref. 9) 

HELM 10 

10 


3.142 

3.144 

20 


6.285 

6.308 

01 


6.285 

6.308 

11 

11 

7.027 

7.027 

12 

12 

12.958 

13.201 

21 

21 

8.889 

8.993 


Table 2. Cutoff Wave Numbers for Circular Waveguide 


Mode 

k c r 

Analytical (ref. 9) 

HELM 10 

TE 01 

3.832 

3.858 

TEn 

1.841 

1.843 

te 12 

5.331 

5.379 

TE 2 i 

3.054 

3.068 

TE 22 

6.706 

6.688 

TM 01 

2.405 

2.413 

TM n 

3.832 

3.854 

tm 12 

7.016 

7.136 

TM 2 i 

5.135 

5.210 

TM 22 

8.417 

8.668 
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Figure 7 . Electric field distribution of some modes for circular waveguide. 


Table 3. Cutoff Wave Numbers for Coaxial Line With T 2 /t\ = 4 


Mode 

k c r i 

Analytical (ref. 10) 

HELM10 

TE n 

0.411 

0.412 

te 21 

0.752 

0.754 

TE 3 i 

1.048 

1.055 

TMoi 

1.024 

1.030 

TMn 

1.112 

1.120 



Figure 8. Cross section of coaxial line. 

2.1.5. Summary 

In section 2, a nodal- based two-dimensional finite 
element method has been described for homogeneous 
waveguides using Galerkin’s technique. The proce- 
dure outlined here is valid for any arbitrary cross 
section of the waveguide filled with homogeneous ma- 
terials. The HELM10 computer program gives the 
cutoff w r ave numbers and the transverse electric field 
components for any mode of propagation in such a 
waveguide. Examples of a rectangular waveguide, a 
circular waveguide, and a coaxial line are presented 
to validate the computer code. Accuracy of the nu- 
merical results depends on the number of elements 
used to represent the geometry. 

2.2. Inhomogeneous Waveguides — Vector 

Formulation 

For inhomogeneous waveguide problems, the 
scalar potential approach is not applicable if one is 


interested in the dispersion characteristics of the 
waveguide. For this a vector finite element approach 
has to be followed. Unfortunately, the nodal-based 
finite element method when applied to the vector 
wave equation results in nonphysical or spurious 
solutions, which are generally attributed to the lack 
of enforcement of divergence condition. (See refs. 2 
and 11.) Many attempts have been made to avoid 
these spurious modes by different variations of the 
nodal-based finite element method. Fortunately, a 
revolutionary approach has been discovered recently. 
This approach uses the so-called “vector basis” or 
“vector elements,” which assign degrees of freedom 
to the edges rather than to the nodes of the element. 
For this reason they are popularly referred to as 
“edge elements.” 

Although the edge elements were first described 
by Whitney (ref. 12) 35 years ago, their importance 
in electromagnetics was not realized until the late 
1980’s (ref. 3). In section 2.2, we describe the for- 
mulation of the finite element method for vector 
wave equation using edge elements. This formula- 
tion, in general, follows the one given by Lee, Sun, 
and Cendes (ref. 13). But instead of attempting to 
solve the dispersion problem right away, we use the 
following steps to build the edge element formulation: 

1. Solve for the homogeneous waveguide prob- 
lem by using two- component transverse vector 
fields 

2. Solve for the generalized waveguide problem 
by using three-component vector fields 

3. Solve for the wave number k Q when the prop- 
agation constant j3 is specified 

4. Solve for the propagation constant /3 of any 
generalized waveguide for a given frequency of 
operation 
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(a) TEn mode. 


(b) TE21 mode. 




(c) TE31 mode. 


(d) TM01 mode. 



(e) TMn mode. 



(f) TM21 mode. 


Figure 9 . Electric field distribution of some modes for coaxial line. 



2.2. 1. Solution of Homogeneous Waveguide 
Problem With Two- Component Transverse 
Vector Fields 

For a homogeneous waveguide, the propagating 
modes can be divided into transverse electric (TE) 
or transverse magnetic (TM) modes, which can be 
solved separately. For the TE mode E z = 0, the 
transverse electric field vector E* satisfies the vector 
wave equation 


v< x (^r V{ x Et ) _ k c £rEt = 0 (42) 

and for the TM mode H z = 0, the transverse electric 
field vector satisfies the vector wave equation 

X (~Vt x H t ) - k 2 c (i r Ht = 0 (43) 


equation (44) can be written in weak form as 


J j — (Vt X Tt) • (Vi X E t)ds = k 2 c e r J I TfEtds 
v r 

- f Tf(nx —Vt x E t ) dl (48) 

dr 7 


On a perfect electric conducting (PEC) boundary, 
the contour integral vanishes as Tt is set to zero to 
satisfy the Dirichlet boundary conditions. Hence, for 
all the problems closed by PEC walls, the last term 
on the right-hand side of equation (48) is set to zero. 
Thus, the final equation can be written as 


j f -^(Vf X Tt) • (Vf x Et) ds = k 2 e r J J T t 
r r 


ds 

(49) 


where p r and E r are the permeability and permittiv- 
ity, respectively, of the material in the waveguide. 

In this section, we illustrate the procedure for the 
TE mode and the same procedure can be applied 
for the TM mode, except the Neumann boundary 
conditions are applied at PEC boundaries instead of 
the Dirichlet boundary conditions. 

2. 2. 1.1. Formulation. For the TE mode, recon- 
sider the wave equation (42), which is 

V< x (Ty t x E() - k 2 £ r Et = 0 

Dot multiplying equation (42) by a vector testing 
function Tt and integrating over the cross section of 
the waveguide, we get 

Jf (r t ■ Vt x (i-V X Et) - k 2 e r Tt • Et) ds = 0 

(44) 

From the vector identities, 

T t • (Vt x A) = (Vt x Tt) • A - Vt • (T t x A) (45) 

(Tt x A) • n = -T t (fix A) (46) 

and the divergence theorem, 

Vt • (Tj x A) ds = j(T t x A )-ndl (47) 

r dr 


2. 2. 1.2. Discretization. The finite elements used 
in section 2.1 are scalar and have unknown parame- 
ters, the values of the scalar field at the nodes of the 
element. Nodal-based finite elements are not suitable 
to represent vector fields in electromagnetics, as the 
boundary conditions often take the form of a specifi- 
cation of only the part of the vector field that is tan- 
gent to the boundary. With nodal-based elements, 
the physical constraint must be transformed into lin- 
ear relationships between the Cartesian components, 
and at nodes where the boundary changes direction, 
an average tangential direction must be determined 
first. These are very difficult (if not impossible) to 
implement in nodal- based finite elements. The fail- 
ure to implement proper conditions results in spuri- 
ous modes which are nonphysical. 

The most elegant and simple approach to elimi- 
nate the disadvantages of the nodal-based elements 
is to use edge elements. Edge elements are recently 
developed finite element bases for vector fields. With 
edge elements, only the tangential continuity of the 
vector fields is imposed across the element bound- 
aries. The advantages of edge elements are as follows: 

1. The edge elements impose the continuity of 
only the tangential components of the electric 
and magnetic fields, which is consistent with 
the physical constraints on these fields. 

2. The interelement boundary conditions are 
automatically obtained through the natural 
boundary conditions. 

3. The Dirichlet boundary condition can be eas- 
ily imposed along the element edges. 
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4. As the edge elements are chosen to be diver- 
gence free, the spurious nonphysical solutions 
are completely eliminated. 

For a single triangular element shown in figure 10, 
the transverse electric field can be expressed as a 
superposition of edge elements. The edge elements 
permit a constant tangential component of the basis 
function along one triangular edge wiiile simultane- 
ously allowing a zero tangential component along the 
other two edges. Three such basis functions over- 
lapping each triangular cell, provide the complete 
expansion (ref. 11): 

3 

E i = (50) 

m=l 

where 

w tm = Lt m (otiVtOij - OcjVtati) (51) 

oti is the first-order shape function associated with 
nodes 1, 2, and 3 defined by equation (15); and Lt m 
is the length of edge m connecting nodes i and j, that 
is, explicitly, the basis functions representing edges 1, 
2, and 3 with coefficients e*i, and e ts written as 

W tl = Lti(a\V tot2 - a^V^ai) (52) 

W*2 = T*2(a2V*a3 — c*3 v i a 2) (53) 

W^3 = L t ^(asVta\ - aiV^) (54) 

The three unknown parameters are etm- It can be 
shown that 

ttm * E i = e tm (on edge m) (55) 

where ttm is a unit vector along edge m in the 
direction of the edge; for example, 

£ ti _ (X2 - Xi)x + (y2 - yi)y 


for edge 1 connecting nodes 1 and 2. In other words, 
etm controls the tangential field on edge m . Later 
Vf • W tm = 0 is verified and hence the electric field 
obtained through equation (50) satisfies V* • E$ = 0, 
the divergence equation, within the element. There- 
fore, the finite element solution is free of spurious 
solutions. 

With the use of the simplex coordinates as de- 
fined by equation (15), the basis function given by 
equation (51) can be written as 

W tm = + B m y)x -f (C m + D m x)y \ (56) 

where 

Am = o>ibj &jbi (57) 

Bjji = fybj Cjbj ( 58 ) 

Cm — &lCj djCi ( 59 ) 

— b t Cj bjC z = B m (60) 

From equation (56), * W* m = 0, which results in 

a divergenceless electric field; that is, V* • E* = 0. 

2. 2. 1.3. Finite element formulation. Substituting 
equation (50) into equation (49) for a single triangu- 
lar element gives the following equation: 

3 

y ^ (V t ^ Win) ’ {&t X 'Wtn)€tm ds 
m~l 

= k 2 c e r If Y. ( W tm ' W tn )etm ds 

A m=l 

(n = 1, 2, 3) (61) 

where A indicates the integration over the triangular 
element. By interchanging the integration and sum- 
mations, equation (61) can be written in a matrix 
form as 

[S el ][e f ] = *2[T el ][e t ] (62) 

Hence the finite element matrices for a single element 
are given by 

[Set] = —J /( V < x w ^) • ( V t x W tn) ds (63) 
A 

and 

[T el ] = e r J j( W tm • W tn ) ds (64) 

A 
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These element matrices can be assembled over all 
the triangles in the cross section of the waveguide to 
obtain a global eigenmatrix equation. 

[S]N = fc c 2 [T][e t ] (65) 

2. 2. 1.4- Finite element matrices . Closed-form 
expressions for integrals in equations (63) and (64) 
can be written as 

//(V t x Wtm) ■ (Vt X W tn ) ds = t^H DmDn 
A 

( 66 ) 

and 

J I (W tm -W tn )ds = ^^(Ia+I t2 +It3+lH+lt5) (67) 

A 

where 

hi = |-( A m A n + C m Cn) J j dx dy (68) 

A 

h2 = ^{CrriDri + C n D m ) j J x dx dy ( 69 ) 

A 

hs = 4~ A n B m ) J J y dx dy (79) 

A 

hi = j(BmBn) j J y 2 dx dy (71) 

A 

% = \{DmD n ) J j x 2 dx dy (72) 

A 

These equations can be reduced further by the 
following integration formulas for integrating over a 
triangle given in reference 14: 


hi — ( A m A n + CmCn) 

(73) 

h2 — {CmF)n 4“ C n D m )x t r i 

(74) 

7*3 = {AmBn + A n Bm)yi T i 

(75) 


(76) 

/« = " (i>, 2 + 94 ,) 

(77) 


2.2. 1.5. Numerical examples. A computer pro- 
gram HELMVEC was written to solve the eigenvalue 
problem presented in section 2.2.1. Numerical data 
computed with this program for a rectangular wave- 
guide and a circular waveguide are given as follows. 

Rectangular waveguide: The cutoff wave numbers 
of an air-filled rectangular waveguide (shown in fig. 4, 
a r /b T = 2) were calculated with the procedure in sec- 
tion 2.2.1 (HELMVEC) and are presented in table 4 
for both the TE and TM modes along with analyt- 
ical data from reference 9. Four hundred triangular 
elements have been used to represent the geometry. 


Table 4. Cutoff Wave Numbers for Rectangular Waveguide 




Mode 

Analytical (ref. 9) 

HELMVEC 

TE mode 

TE 10 

3.142 

3.140 

TE20 

6.285 

6.274 

TE 30 

9.428 

9.396 

TE 01 

6.285 

6.274 

TE„ 

7.027 

7.024 

te 21 

8.889 

8.897 

TE 31 

11.331 

11.339 

TEq2 

12.571 

12.497 

TM mode 

TMn 

7.027 

7.024 

TMl2 

8.889 

8.889 

TM 13 

11.331 

11.340 

TM 21 

12.958 

12.902 

TM 22 

14.055 

14.044 

TM 23 

15.710 

15.757 

TM 31 

19.117 

18.906 

TM 32 

19.877 

19.728 


Circular waveguide: The cutoff wave numbers for 
an air-filled circular waveguide (Radius = 1) were 
also computed by HELMVEC and compared with 
analytically calculated values from reference 9. These 
results are presented in table 5 for both the TE and 
TM modes. The waveguide geometry is represented 
by 200 triangular elements. 

2.2.2. Inhomogeneous Waveguide Problems 

Using Three- Component Vector Fields 

In this section, a generalized approach for find- 
ing the cutoff frequencies of inhomogeneously filled 
waveguides is described. 
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Table 5. Cutoff Wave Numbers for Circular Waveguide 



k c 

Mode 

Analytical (ref. 9) 

HELMVEC 


TE mode 


TEoi 

3.832 

3.834 

TE U 

1.841 

1.846 

te 12 

3.054 

3.070 

TEi 3 

4.201 

4.192 

TE 20 

7.016 

7.002 

te 21 

5.331 

5.332 

te 22 

6.706 

6.620 

TE 23 

8.015 

8.178 


TM mode 


TM 10 

2.405 

2.416 

TM U 

3.832 

3.831 

tm 12 

5.136 

5.155 

tm 13 

6.380 

6.358 

TM 20 

5.520 

5.508 

TM 21 

7.016 

7.054 

tm 22 

8.417 

8.403 

tm 23 

9.761 

9.644 


and equation (81) can be written in its weak form by 
following equations (2) through (7) to obtain 

J / ' VtEz ^ ds = ^ / J TzEz ds ( 83 ) 

r r r 

where the testing function is T = Tt + zT z . 

2. 2. 2. 2. Discretization. As the vector Helmholtz 
equation is divided into two parts, vector-based tan- 
gential edge elements can be used to approximate 
the transverse fields, and nodal-based first-order 
Lagrangian interpolation functions can be used to 
approximate the z-component. FYom equation (50), 
which is 

3 

~ y ^ w i m 

m= 1 

where m indicates the rath edge of the triangle and 
W tm is the edge element for edge ra. The testing 
function 7J is chosen to be the same as the basis 
function in equation (50) (shown above); that is, 
Tt — W t m . The z-component of the electric field 
can be written as 


2.2.2. 1. Formulation. This generalized approach 
can be followed by using either the E or H field. We 
will illustrate the case for the E field. The vector 
wave equation for E is given by 

v< X x E ) _ k c £ r E = 0 (78) 

The electric field can be written as 

E = E t + zE z (79) 


3 

= £ e zi Q i (84) 

*=1 

Here i indicates ith node and a* is the simplex 
coordinate of node i as given in equation (15). Also 
the testing function T z is chosen to be the same as 
the basis function in equation (84); that is, T z = aj. 

2.2. 2. 3. Finite element formulation. Substituting 
equations (50) and (84) into equations (82) and (83), 
respectively, we can get the following equations: 


Hence equation (78) can be divided into two parts — 
one consisting of the transverse electric fields and the 
other, the z-component of the electric field — that is, 

Vt x (^ v< x e 0 - k2£rEt = ° (8o) 
v t ■ —V t E z + k 2 e r E z = 0 (81) 

fl r 

Equation (80) can be written in its weak form as 

/ / ^ (v< x Tt) ■ x Et)ds = *&• / 7 Tt • Etds 

r r 

(82) 


3 

Me (Vt x W tm ) • (Vi x W tn)etmds 
r A m=1 

= k c £ r j J ) ^ • W tn )e t mds 


A m=1 
(n= 1,2,3) 


(85) 


^£ e zi J f ^ t0ii '^ t0t ^ dxd y 

t=i A 

^ €zi 1 1 ai< * J dX dy ^ = 2 ’ ^ (®6) 


— k c e T 


13 




By interchanging the summation and integration in 
equation (85), these two equations can be combined 
to be written in a matrix form as 


1 

Co 

CD 

o 

1 


CN < 
-Sc 

II 

T el(t) 

1 

o 

’e t m 

t 

o 

s 

1 

Le z J 


o 

T el (*) J 



Circular waveguide: The cutoff wave numbers for 
an air- filled circular waveguide (Radius = 1) are also 
computed and compared with analytically calculated 
values. These are given in table 7. 

Table 6 . Cutoff Wave Numbers for Rectangular Waveguide 


Hence the element matrices are given by 
Sei W = ~ r jj(V t X Wftn) • (V* X Wtn) ds (88) 


^el(z) J J ) dx dy 

A 

■^el(i) ~ ~ r f J ds 

A 

= £ r II otictj dx dy 


l el(z) 


(89) 

(90) 

(91) 


These element matrices can be assembled over all 
the triangles in the cross section of the waveguide to 
obtain a global eigenmatrix equation: 


- St 

0 ■ 

m et* 

= k 2 

' T t 

0 ‘ 

'e t ~ 

. 0 

Sz. 



. 0 

T z . 

. £ z . 


(92) 


Mode 


Analytical (ref. 9) 

HELMVECl for— 

E 

H 

TE mode 

TEio 

3.142 

3.142 

3.144 

TE 20 

6.285 

6.274 

6.308 

TE 30 

9.428 

9.396 

9.511 

TE 01 

6.285 

6.274 

6.308 

TE U 

7.027 

7.024 

7.024 

TE 2 i 

8.889 

8.897 

8.897 

TE 31 

11.331 

11.339 

11.340 

TEq 2 

12.571 

12.497 

12.772 

TM mode 

TM„ 

7.027 

7.072 

7.207 

TM 2 1 

8.889 

8.995 

8.993 

TM 31 

11.331 

11.541 

11.537 

TMq 2 

12.571 

12.498 

12.772 


Table 7. Cutoff Wave Numbers for Circular Waveguide 


2. 2. 2-4- Finite element matrices. The finite ele- 
ment matrices shown in equations (88) through (91) 
are already derived in the previous sections. Equa- 
tion (88) is given by equation (66), equation (89) by 
equation (31), equation (90) by equation (67), and 
equation (91) by equation (33). 

2.2.2. 5. Numerical examples. A computer pro- 
gram HELMVEC1 was written for the three- 
component vector formulation to calculate cutoff 
wave numbers of waveguides with arbitrary cross sec- 
tions. Numerical data computed with HELMVECl 
for a rectangular waveguide and a circular waveguide 
are given as follows. 

Rectangular waveguide: The cutoff wave numbers 
of an air- filled rectangular waveguide (shown in fig. 4, 
a r /b r = 2) are calculated with the procedure in sec- 
tion 2.2.2 (HELMVECl) and are presented in table 6 
along with analytical data from reference 9. In the 
present formulation, it is not required to calculate 
eigenvalues separately for TE and TM modes. Both 
E and H formulations result in almost identical nu- 
merical results. 


Mode 

k c 

Analytical (ref. 9) 

HELMVECl for— 

E 

H 

TE mode 

TE 01 

3.832 

3.798 

3.891 

TE„ 

1.841 

1.869 

1.875 

TE 12 

3.054 

3.092 

3.172 

te 13 

4.201 

4.258 

4.084 

TE 20 

7.016 

6.899 

6.912 

te 21 

5.331 

5.426 

5.226 

TE 22 

6.706 

6.653 

6.746 

TE 23 

8.015 

7.930 

8.116 

TM mode 

TM 10 

2.405 

2.481 

2.439 

TMji 

3.832 

3.798 

3.891 

TMl 2 

5.136 

5.228 

5.183 

TM 13 

6.380 

6.507 

6.335 

TM 20 

5.520 

5.493 

5.551 

TM 21 

7.016 

7.490 

7.410 

TM 22 

8.417 

8.591 

8.449 

TM 23 

9.761 

9.769 

9.804 
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2.2,3. Wave-Number Determination for Given 
Propagation Constant 

A generalized approach for finding the wave 
number k 0 is described for a given propagation 
constant p. 


~JJ T ^t ■ &tE z + E t )] ds = k 2 0 e r J f T Z E Z ds 
r r 

( 101 ) 

With the vector identities 
A • (V* x B) = (V t x A) ■ B - V t ■ (A x B) (102) 


2.2.3. 1. Formulation. Again the same approach 
as used earlier can be followed by using either the 
E or H field. We illustrate the case for the E field, 
which is the same for the H field. The vector wave 
equation for the E field is given by 

V X (~V x e) - fc 2 e r E = 0 (93) 

where E = ( E x x + E y y 4- E z z) exp(-jfiz). By doing 
the curl-curl operation and separating the transverse 
from the longitudinal components, equation (93) can 
be divided into two equations and rewritten as 

Vt X (~ V * X Et J - — (; jdV t E z - /3 2 E = k 2 0 e r E t 

} (94) 

— -[Vt • (V t E z + j/3E f )] = k 2 0 e r E z (95) 
fi T 

and Et = E x x + E y y. Separate 0 from and to 
have real- valued matrices introduce the scaling 

E t 

Et = j (96) 

E z = ~Y ( 97 ) 

Then equations (94) and (95) can be written as 


V( X ( — V t xE,) + i (0 2 V t E z + P 2 E,) = kle r E t (98) 

\Pr J Mr 


[V t • (V t E z + Et)] = k 2 e r E z (99) 

fl r 

To apply Galerkin’s method to equations (98) 
and (99), multiply equation (98) with the testing 
function T t and equation (99) with the testing func- 
tion T z and integrate both the equations over the 
cross section of the waveguide T; that is, 


Hfr 


( — V( X Et) + —(Tt ■ V t E z + Tf Et 

\Pr ) fir 


ds 


= k 2 0 e r 



* E t ds 


( 100 ) 


/ f Vt • (A x B) ds = J(AxB)ndl = - J A(nxB )dl 
r dr dr 

(103) 

V t - /A = A ■ V</ + /V t • A (104) 

//v,.a*-/a •ndZ (105) 
r dr 

equations (100) and (101) can be written in their 
weak form as 

^ II f (Vt x r<) ' (v ‘ x Ei) + O 9 * 3 * ' VEz + ^ Tt ' E *)] ds 

v 

= k o £ r I j TfE t ds-j- j Tf(nxV xEt)ds (106) 


Pff WtT z ■ V t E z + V t T z • Et) ds 

J 

r dr 

(t z ~ + T z n ■ Et^j ds (107) 


— fcp£r J j T Z E Z ds + — 


If the waveguide boundary dT is assumed to be 
perfectly conducting, then 7} = 0 and T z = 0 on 
dT. Hence, the line integrals on the right-hand side 
of equations (106) and (107) can be neglected. Mul- 
tiplying equation (107) with /? 2 for the sake of sym- 
metry, equations (106) and (107) can be rewritten 
as 


i // [<v, x T,MV, X E,) 
r 

+ [p 2 TfVE z + p 2 T t -E t 'j]ds 

= k 2 0 e r j j TfE t ds (108) 

r 
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(109) 


— If (V t T z ■ V t E z + V t T z ■ E t ) ds 
T 


T Z E Z ds 


2. 2.3.2. Discretization. As the vector Helmholtz 
equation is divided into two parts, vector-based tan- 
gential edge elements can be used to approximate 
the transverse fields, and nodal-based first-order 
Lagrangian interpolation functions can be used to 
approximate the z-component. From equation (50), 
which is 

3 

Et ~ ^ tm 

m— 1 

where m indicates the mth edge of the triangle and 
W^ m is the edge element for edge m. The testing 
function Tt is chosen to be the same as the basis 
function in equation (50); that is, Tt = W* m , The 
z-component can be written as 


Ej z — ^ ^ ^zi^i 

1 

(This equation is eq. (84).) Here i indicates zth node 
and ai is the simplex coordinate of node i as given in 
equation (15). Also the testing function T z is chosen 
to be the same as the basis function in equation (85) 
(given above); that is, T z = a*. 

2.2.3. 3. Finite element formulation. Substitut- 
ing equations (50) and (84) into equations (108) 
and (109), respectively, integrating over a single tri- 
angular element, and interchanging the integration 
and summation give 

3 

— E //(V f x W fm ) • (Vf x W tn )e tm ds 
iV 

+ — E l f ( W tm ‘ Vaj)e zj ds 

Mr -= lV 

+ — E [ W tn )e tm ds 

MT m=li J 

— kl 'y ; £r j J '(Wtm ' W<n) e tm ds 


m=l A 

(n = 1,2,3; j = 1,2,3) 


L 

Mr 


E [ f ( Va » ' ds 

i = 1 A 

g2 3 . . 

H E // (^ Q * ' ds 

i=l V 

= kl f3 2 £ T j j aiOcje z i ds 


(j= 1,2,3; n = 1,2,3) 


( 111 ) 


Subscripts for q indicate node numbers and sub- 
scripts for Wf indicate edge numbers. Equa- 
tions (110) and (111) can be written in matrix form 
as 


1" ^el (tt) 

Co 

1 

"e* " 

= ki 

^el(tt) 0 1 

~et~ 

-K> 

CO 

i 

^ e\(zz ) J 

.e z \ 


o 

£3 

TT 

tt 

i — 

.e z . 


The element matrices are given by 

S em = — r IJ(V t x W tm) ■ (Vt X W tn ) ds 


( 112 ) 


+ — r I f (W tm . W tn ) ds 
A 

(113) 

{ U) = J T IJ^-^ds 

A 

(114) 

S el(zt) = ^;Jf(Va i -W tn )ds 

A 

(115) 

5 el( Z2 ) = Jj ' Va >) 

r A 

(116) 

^el(tt) ~ Sr j J ^ tn ) 

A 

(117) 

^el(zz) = P^ e r J J OLiOtjds 

(118) 


(110) 


These element matrices can be assembled over all 
the triangles in the cross section of the waveguide to 
obtain a global eigenvalue equation as follows: 


(119) 


-s tt 

Stz 

m et' 

= ki 

■ Ttt 

0 ' 

' et~ 

.s zt 

Szz. 

.e z . 

0 

. 0 

T zz . 

e z . 
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2. 2. 3. 4- Finite element matrices . Closed-form 
expressions are derived for the finite element matri- 
ces in equations (113) through (117). From equa- 
tions (66) and (67), 


5el ^> ~ 


+ X> 

V fc=i t 


(120) 


With equations (73) through (77) and the proce- 
dures in sections 2.2.1 and 2.2.2, the remaining ele- 
ment matrices are derived as follows: 


02 

^el(tz) = ~ ^2 i^m + ^m2/tri)+ c j{Cm + 7)m^tri)] 


fl r 

( 121 ) 

02 

^el(zt) ^ g^2 B n y^ r j)-|- Ci(C n + -Dri^tri)] 


0 1 bibj + CiCj 

S e\ (zz) = 


fir 4A 


rji ^ FtmLt n j 

I e\{tt) - e r lgA 3 J tk 


T e\ (zz) = /^r — 


= /3 Z S r 


2. A 
” 12 


(m = n) 
(m / n) 


(122) 

(123) 

(124) 

(125) 


2.2.3. 5. Numerical example. A computer pro- 
gram HELMVEC2 was written to calculate the wave 
number of any generalized waveguide structure for a 
given propagation constant. To check the validity of 
the procedure, the first 10 wave numbers with /5 = 10 
for LSM modes of a partially filled square waveguide 
(fig. 11) have been obtained by HELMVEC2 and 
compared with the available results in the literature 
(ref. 15). (See table 8.) The geometry is modelled by 
72 triangular elements. 
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Figure 11. Partially filled square waveguide. 


Table 8. Wave Numbers for LSM Modes of Square 
Waveguide With 0 = 10 


Mode 

hoL 

HELMVEC2 

Hayata et al. (ref. 15) 

1 

8.815 

8.8093 

2 

9.443 

9.3896 

3 

10.35 

10.2752 

4 

11.141 

11.103 

5 

11.2890 

11.2677 

6 

11.4246 

11.4501 

7 

12.146 

11.9882 

8 

12.5894 

12.6686 

9 

12.8237 

12.8092 

10 

12.9987 

12.9575 


2.2.4 • Dispersion Characteristics of Waveguides 

This section is simply an extension of sec- 
tion 2.2.3. The finite element equations are re- 
arranged to obtain j3 when the operating frequency 
(or the wave number) is specified. 

2.2.4- 1- Formulation. The weak form derived in 
section 2.2.3. 1 can be rearranged and equations (108) 
and (109) can be written as 

/v //( Vf x ' ^ Vt x Et ^ ds ~ k ° £r JJ TtEtds 

r r 



= k 2 0 f3 2 e r j f T z E z ds (127) 

r 


2.2.4 -2. Discretization. As the vector Helmholtz 
equation is divided into two parts, vector-based tan- 
gential edge elements can be used to approximate 
the transverse fields, and nodal-based first-order 
Lagrangian interpolation functions can be used to 
approximate the z-component. From equation (50), 
which is 

3 

m= 1 
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where m indicates the mth edge of the triangle and 
W t m is the edge element for edge m. The testing 
function Tt is chosen to be the same as the basis 
function in equation (127); that is, T± = Wf m . The 
z-component is written as 


E z = 


3 


z / v e z {Or % 
i= 1 

(which is eq. (84)). Here i indicates ith node and 
a* is the simplex coordinate of node i as given in 
equation (15). Also the testing function T z is chosen 
to be the same as the basis function in equation (50) 
(given above); that is, T z — a*. 

2.2. 4-3. Finite element formulation. Substitut- 
ing equations (127) and (50) into equations (125) 
and (126), respectively, integrating over a single tri- 
angular element, and interchanging the integration 
and summation give 


3 

^ E / / ( v * x W tm ) - (Vt X W tn )etm 
T 771=1 A 

xds-kl^er / / (W t m ’ ds 

m= 1 V 

]_ 3 ^ r r 

— 'y ^ J J (AV^ m • ^7ckj)e z j ds 
r m=1 A 

^E //(W im -W tn )e tm ds 


= - 0 * 


(n — 1, 2, 3; j — 1, 2, 3) 


(128) 


P_ 

Mr 


|2 3 ^ » - 

- E J j ( Va * ' Va j) e zi 

t=l A 


+ ^E ff(V<* i -'Wtn)etnda 
r 7=1 A 

ft ^ ^ J' j &iOtj€zi ds 


(7 = 1,2,3; 71=1,2,3) 


(129) 


Subscripts for a indicate node numbers and sub- 
scripts for W t indicate edge numbers. Equations (84) 


and (128) can be written in matrix form as 


■^el(tt) -^el (tz) 


o 

CO 

'e t ‘ 

= -J 

0 OJ 

e z . 



l T e\{zt) T e\(zz)\ 


et 
e z J 
(130) 


The element matrices are given by 

^e\(tt) = ~ ff(Vt x w tm ) • (Vt X W tn ) ds 
A 

-k 2 0 e r ff(W t m-Wtn)ds ( 131 ) 

A 

T e \(tt) = £ r ff(Wtm-Wtn)ds ( 132 ) 

A 

T<A(tz) — — J J (Wtm ■ VQj) ds ( 133 ) 

r A 

^el(^) = E jj (Vat • Wtn) ds ( 134 ) 

r A 

T e\(zz) = J j f ( Va * ' Vo?) ds - kle r J j aiCkj ds 
A A 

( 135 ) 

These element matrices can be assembled over all 
the triangles in the cross section of the waveguide to 
obtain a global eigenvalue equation. 


et 


Stt 0 

0 OjLe, 


= (~P 2 ) 


Ttt Ttz 

L Fzt T zz 


et 

e z 


(136) 


2. 2.4-4- Finite element matrices. Closed-form 
expressions are derived for the finite element matri- 
ces in equations (130) through (134). From equa- 
tions (66) and (67), 


_ J 


LtmLtn ( ^ ^ 7 2 j 

16 A 3 [ ° mDn " k ° Er ^ Itk 


(137) 

With equations (73) through (77) and the proce- 
dures in sections 2.2.1 and 2.2.2, the remaining ele- 
ment matrices are derived as follows: 

-^el (tz) = ~ g^2 [bj(Am + B m y iT i) 


T T -^m^tri)] 

T e \(zt) = “^2 + BnVtri) 

+ Ci(C n + D n x tr i] 


(138) 


(139) 
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c 1 bfij + <H c j 

el < zz > ~~iT r 4A 


(140) 


rr-i LtmLtn t 

■'el («) ~ £ r 16A 3 Z, 

k=\ 


(141) 


-^el(zz) 


i »? + ■? 

/i r 4A 



1 b-i bj “}~ Ci Cj 

p r 4A 



(*=j) 

> 

(Mj)^ 

(142) 


2.2. 4-5. Numerical examples. A computer pro- 
gram HELMVEC3 was written to calculate the 
dispersion characteristics of inhomogeneously filled 
waveguides. For the first example for a partially filled 
rectangular waveguide with b r /a r = 0.45, d/b r = 0.5, 
and £ t = 2.45 shown in figure 12, the values of f3/k 0 
are calculated and compared with analytical results 
given by Harrington (ref. 9, p. 161). The numeri- 
cal results are presented in table 9 as a function of 
a r j \ 0 . The rectangular waveguide for this first ex- 
ample is modelled by 100 triangular elements. 


Table 9. Dispersion Characteristics of Partially 
Filled Rectangular Waveguide of Figure 12 


brf A? 

0/ko 

Analytical (ref. 9) 

HELMVEC3 

0.2 

0.48 

0.47 

0.3 

1.00 

1.01 

0.4 

1.18 

1.17 

0.5 

1.26 

1.28 

0.6 

1.30 

1.35 



Figure 12. Partially filled rectangular waveguide with 
b T ja r - 0.45, d/br = 0.5, and = 2.45. 


For the second example for a partially filled rect- 
angular waveguide with b T /a T = 0.45 and e r — 2.45 
shown in figure 13, the dispersion characteristics were 
also calculated and compared with those given by 
Harrington (ref. 9, p. 162). The numerical results are 
presented in table 10. The rectangular waveguide for 
this second example is also modelled by 100 triangu- 
lar elements. 




Y 

A 

► X 


Figure 13. Partially filled waveguide with b r /a r = 0.45 and 

5 r = 2.45. 

2.2.5 . Summary 

In section 2.2, a detailed formulation for vector 
finite elements for two-dimensional eigenvalue prob- 
lems in electromagnetics has been presented. Edge 
elements are used to impose divergence free condition 
on fields and hence the nonphysical modes (ref. 2), 
which plagued vector finite elements for many years 
are avoided. When assembling the global matrices 
from element matrices with edge basis functions, a 
unique global edge direction is defined to ensure field 
continuity across all edges. The basis function is mul- 
tiplied by —1, if the local edge vector does not have 
the same direction as the global edge direction. The 
numerical examples presented show the validity and 
accuracy of the analyses and computer codes, respec- 
tively. Both the formulation and the computer codes 
are valid for any arbitrarily shaped waveguides filled 
with inhomogeneous materials. 

3. Three-Dimensional Problems 

3.1. Eigenvalues of Three-Dimensional 

Cavity — Vector Formulation 

The problem of calculating resonant frequencies 
of three-dimensional cavities has been plagued by 
spurious modes for many years. As mentioned in sec- 
tion 2.2, this problem has been recently overcome by 
using vector-based tangential edge elements (refs. 15 
and 16). In section 3.1, we formulate the Galerkin 
finite element method for three-dimensional cavities 
and present results for various shapes of cavities. 
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Table 10. Dispersion Characteristics of Partially 
Filled Rectangular Waveguide of Figure 13 


3.1.1 . Formulation 


dr / 

0/k o 

Analytical (ref. 9) 

HELMVEC3 

d/a T = 0 

0.4 



0.5 

0.03 

0.04 

0.6 

0.52 

0.56 

0.7 

0.70 

0.71 

0.8 

0.79 

0.78 

0.9 

0.83 

0.83 

1.0 

0.88 

0.87 

d/a r = 0.167 

0.4 



0.5 

0.21 

0.18 

0.6 

0.60 

0.59 

0.7 

0.72 

0.74 

0.8 

0.82 

0.81 

0.9 

0.88 

0.87 

1.0 

0.91 

0.90 

dja T — 0.286 

0.4 



0.5 

0.51 

0.44 

0.6 

0.78 

0.74 

0.7 

0.90 

0.88 

0.8 

0.99 

1.05 

0.9 

1.03 

1.03 

1.0 

1.10 

1.09 

d/a r = 0.375 

0.4 



0.5 

0.68 

0.66 

0.6 

0.91 

0.90 

0.7 

1.05 

1.03 

0.8 

1.13 

1.11 

0.9 

1.20 

1.18 

1.0 

1.25 

1.23 

d/a r = 0.5 

0.4 

0.40 

0.42 

0.5 

0.90 

0.89 

0.6 

1.10 

1.09 

0.7 

1.20 

1.19 

0.8 

1.25 

1.24 

0.9 

1.30 

1.31 

1.0 

1.35 

1.35 

d/a T — 0.6 

0.4 

0.70 

0.67 

0.5 

1.02 

1.03 

0.6 

1.18 

1.19 

0.7 

1.23 

1.27 

0.8 

1.31 

1.33 

0.9 

1.38 

1.37 

1.0 

1.41 

1.40 

d/a T = 0.8 

0.4 

0.90 

0.91 

0.5 

1.18 

1.18 

0.6 

1.29 

1.30 

0.7 

1.38 

1.37 

0.8 

1.41 

1.42 

0.9 

1.43 

1.44 

1.0 

1.44 

1.47 


This formulation can be followed by using either 
the E or H field. We illustrate the E field formula- 
tion. Consider the vector wave equation, 


V x 


^VxEj-/c^ r E = 0 (143) 


where 

E = E x x + E y y + E z z 

Introduce a testing function T = T x x + T y y + T z z. 
To apply Galerkin’s method, multiply equation (143) 
with T and integrate over the volume of the cavity v: 


kle r T ■ Ej dv = 0 
(144) 

A • (V x B) = (V x A) - B - V • (A x B) (145) 
equation (144) can be written as 


///( tvx (^ vxe )- 

V 

With the vector identity, 


■ E ds 

V v 

+ /// v l rx (£ vx ’)]* (14< 


Invoking the divergence theorem 


/// v 


■ Adv = 



A - rids 


(147) 


and 

(A x B) • n = -A • (n x B) (148) 

where v indicates the integration over the volume of 
the cavity, S indicates the integration over the outer 
surface of the cavity, and n is the outward unit vector 
normal to the surface, equation (146) becomes 


/// (vxT) -(i vxE )* , -' fc /// T - B 

V V 

-JI T [* *&**)}*“ ( 


dv 


For a cavity bounded by perfectly conducting 
electric conductor, the field as well as the testing 
function T has to be zero on the outer surface; hence, 
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the last term on the right-hand side of equation (149) 
vanishes. Equation (149) can be written in its final 
form as 


///,Vx T) .(iv, E ) 


dv 


-+JIJ TE dv (150) 


3.1.2 . Discretization 

The volume of the cavity is discretized by us- 
ing first-order tetrahedral elements such as the one 
shown in figure 14. The first-order tetrahedron has 
four nodes and six edges. The six edges are formed 
as shown in table 11. 


l 



Figure 14. First-order tetrahedral element. 


Table 11. Formation of Edges of Tetrahedral Element 


Edge, m 

Node 

i 

j 

1 

1 

2 

2 

2 

3 

3 

1 

3 

4 

1 

4 

5 

2 

4 

6 

3 

4 


The electric field in a single tetrahedral element 
is represented as 


6 

E = £ e m W m (151) 

m = 1 

The six unknown parameters associated with each 
edge are e\, e2> • * •, eg. The total field is obtained by 
evaluating equation (151). 


The vector tangential edge elements W m are 
given by (ref. 3) 

"Wm = ocfj — a^Vck^j) (152) 


In equation (152), m stands for edge number, i and 
j stand for the nodes connecting edge i, L m is the 
length of the edge m, a a and a t j are the simplex 
coordinates associated with nodes i and j. 

The simplex coordinates for the nodes of a tetra- 
hedron element are given by Silvester and Ferrari in 
reference 1 


Vi 

a * 1= V 

(153) 

V 2 

<*t2=y 

(154) 

v 3 

<*t3=y 

(155) 

V 4 

at4= V 

(156) 


where V is the volume of the tetrahedron given by 




1 

Xl 

Vl 

Z\ 


1 

1 

X 2 

V 2 

*2 

V - 

— 

1 




6 

x 3 

V3 

z 3 



1 

X 4 

y\ 

z 4 



1 

X 

y 

z 

Vi = 

1 

1 

*2 

y 2 

z 2 

— 

1 



6 

x 3 

2/3 

z 3 



1 

X 4 

yi 

Z 4 



1 

x x 

y\ 

z l 

V 2 = 

1 

1 

X 

y 

Z 





6 

1 

x 3 

y3 

z 3 



1 

X 4 

2/4 

z 4 



1 

xi 

yi 

Z\ 

V 3 = 

1 

1 

x 2 

V 2 

z 2 

— 

1 



6 

X 

y 

z 



1 

X 4 

2/4 

z 4 



1 

Xl 

y\ 

z \ 

V 4 = 

1 

1 

x 2 

V 2 

z 2 

— 




6 

1 

x 3 

y3 

z 3 



1 

X 

y 

Z 


(157) 


(158) 


(159) 


(160) 


(161) 


For any node (f = 1,2, 3, 4), a t i is given by 


Oiti = 


a ti + bfjX -b cgy -f dgz 
6V 


(162) 
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where ba , c^, and da are appropriate cofactors 
picked from the determinants in V \ , V2, ^3, and V4 
for m — 1,2,3, and 4, respectively. From equa- 
tions (152) and (162), the edge elements are given 
as 

L _ 

+ B xm y -h 


+ (^4ym + Bym x + Cymz)y 

+ (A 

znn + B zm x + C zm y)z ] (163) 

where 

= a tjbti (164) 

Bxm = c tibtj c tj^U (165) 

Cxm — dfobij — dfjba (166) 

Aym — a tj c ti (167) 

Bym — btjCfj bfjCa — B X m (168) 

Gym = dactj dijCa (169) 

= ~ Q'tjd'ti (170) 

Bztti ~ badtj bfjda = C X m (1^1) 

Czm ~ dtj Ctjdti = Gym (172) 


Also W m can be shown to satisfy the condition, 


where the element matrices are given by 

[S d ] = — ///(V x W m ) ■ (V x W B ) dv (176) 

▲ 

[Tel) = £r f fI (Wm ' Wn) dV (17?) 

▲ 

These element matrices can be assembled over all the 
tetrahedral elements in the cavity volume to obtain 
a global eigenmatrix equation, 

[S][e] = fcg[T][e] (178) 

To ensure field continuity across all edges, a 
unique global edge direction is defined (i.e., always 
pointing from the smaller node number to the larger 
node number) so that equation (151) has to be mul- 
tiplied by —1, if the local edge vector does not have 
the same direction as the global edge direction. 

The electric field is zero on the PEC bound- 
aries. It is imposed by taking the coefficients of equa- 
tion (151) as zeros. In other words, the edges on the 
boundary are simply ignored when forming the finite 
element matrices and hence reducing the order of the 
matrices to be solved. 


*-' W — ft (othS) (173) 

where t m is the unit vector along the direction of the 
edge. 


3.1.4 • Finite Element Matrices 

The aim of this section is to obtain closed-form 
expressions for equations (176) and (177). From 
equation (163), 


3.1.3. Finite Element Formulation 


Substituting equation (151) into equation (150), 
integrating over the volume of one tetrahedral ele- 
ment, and interchanging the summation and integra- 
tion give 


1 

fi r 



(V x W m ) • (V x W n )e m dv 


6 

= k l E £r ///(W m ■ W n )e m dv 

771 — 1 a 


V X Wra — Jgj/2 ( Gzm x + C X mV + By m Z ) (179) 

and hence 

(V x W m ) - (V X W„) = ^^2 

X ( C Z mCzn + C X mCxn + By m By n ) (180) 

From equations (176), (177), and (180), and mak- 
ing use of the integration formulas given in refer- 
ence 18, the closed-form expressions for element ma- 
trices are given as 


(n = 1,2, . .. ,6) (174) 

where A indicates integration over the volume of 
tetrahedron. This can be written in matrix form as 

[5 el ][e] = fc 0 2 [r el ][e] (175) 


? e l — 324 V- 3 /L 6 {C zrn C Z n + C X mCxn + Bym By n ) 

(181) 

(182) 


k = 1 
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where 

h — AxmAxn A Ay m Ayn + A zm A zn 

I2 = ( Ay m Byn A AynBym + A zm B zn + A zn B 2?7l ) Xtet 
^3 = (-^xm^xn 4 " - 4 .xn-Bxm 4 " -A^mCzn 4 - ^ 4 271^2771) 2 /tet 

^4 = (-^4x771^X71 4“ AxnCxm 4“ Ay m Cyn A Ay n Cym) Z^Qt 
h = ^{^zmCzn A B zn C zm ) ^iVt 4- 16xtetS/tet^ 

^6 = BxmCxn A B xn C X Tn ) Vi^i 4~ 16t/ tet Ztet^J 

h — ^{BymCyn 4" By n Cy m ) %iZi 4- 16xtet^tet^j 

^8 = 2Q {BymByn A B zm B zn ) x\ A 16^tet^ 

= 2Q {BxmBxn + C zm C zn ) (V Vi A 162/tet^j 

ho = 2q(Cx7t \C xn + CymCyn) zf 4- 16ztet^j 


3.1.5. Numerical Examples 


Z 



Figure 15. Air-filled rectangular cavity. Size: 1 by 0.5 by 
0.75 cm. 


Y 



Figure 16. Half-filled rectangular cavity with dielectric mate- 
rial e r = 2.0 and filled from z = 0.5 to 1 cm. Size: 1 by 
0.1 by 1 cm. 


A computer program FEM3D0 was written to cal- 
culate the eigenvalues of a three-dimensional cavity. 
For a three-dimensional problem, the number of vari- 
ables increases drastically compared with those for a 
two-dimensional problem. Hence it is not economi- 
cal to use a generalized eigenvalue solver. The pro- 
gram FEM3D1 was written to take advantage of the 
sparse nature of the finite element matrices. This 
program exploits the symmetry of the matrices and 
stores only the nonzero entries in the lower triangle 
of the matrices; hence, considerable savings in mem- 
ory result. It also makes use of the sparse eigenvalue 
solvers available in VECLIB (ref. 19) and results in 
faster computation. 

The following numerical examples have been 
taken from Chatterjee, Jin, and Volakis (ref. 17), in 
which the resonant wave numbers for various three- 
dimensional cavities have been calculated by using 
a different kind of edge elements. The analytical re- 
sults mentioned here are also taken from reference 17. 

Air-filled rectangular cavity: The eigenvalues of 
an air-filled rectangular cavity (shown in fig. 15) 


were calculated by FEM3D1 and are presented in 
table 12. The cavity geometry is represented by 
343 tetrahedral elements. 


Table 12. Eigenvalues of Air-Filled Rectangular Cavity 


Mode 

k Q} cm 1 

Analytical 
(ref. 17) 

FEM3D1 
(343 elements) 

Reference 17 

TE 101 

5.236 

5.242 

5.213 

TMno 

7.025 

6.942 

6.977 

TEon 

7.531 

7.372 

7.474 

TE 201 

7.531 

7.560 

7.573 

TE m 

8.179 

8.064 

7.991 

TM m 

8.179 

8.164 

8.122 

TM 210 

8.886 

8.725 

8.572 

TE 102 

8.947 

8.871 

8.795 


Half-filled rectangular cavity: A half-filled (filled 
from z = 0.5 to 1 cm) rectangular cavity with 
dielectric material e r ~ 2.0 is shown in figure 16. 
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The eigenvalues were calculated by FEM3D1 and 
are presented in table 13. The cavity geometry was 
modelled with 615 tetrahedral elements. 


Table 13. Eigenvalues of Half-Filled Rectangular Cavity 


Mode 

k Q , cm 1 

Analytical 

FEM3D1 
(615 elements) 

Reference 17 

TE^ioi 

3.538 

3.524 

3.534 

TE Z 201 

5.445 

5.401 

5.440 

TE.I02 

5.935 

5.931 

5.916 

TEz301 

7.503 

7.382 

7.501 

TE z 202 

7.633 

7.562 

7.560 

TE z 103 

8.096 

8.003 

8.056 


Circular cylindrical cavity: The eigenvalues of an 
air-filled circular cylindrical cavity (shown in fig. 17) 
were calculated by FEM3D1 and are presented in 
table 14. The geometry of this cavity was represented 
by 633 tetrahedral elements. 



Figure 17. Air-filled circular cylindrical cavity. Dimensions 
are in centimeters. 


Table 14. Eigenvalues of Air-Filled Circular 
Cylindrical Cavity 


Mode 

k Ql cm 1 

Analytical 

FEM3D1 
(633 elements) 

Reference 17 

TMoio 

4.810 

4.782 

4.809 

TEni 

7.283 

7.210 

7.202 

TE m 

7.283 

7.290 

7.288 

TMno 

7.650 

7.575 

7.633 

TMno 

7.650 

7.620 

7.724 

TMon 

7.840 

7.901 

7.940 

TE 2 n 

8.658 

8.676 

8.697 

TE 211 

8.658 


8.865 


Spherical cavity: The spherical cavity geometry 
was represented by 473 tetrahedral elements. The 
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eigenvalues, calculated by FEM3D1, are presented in 
table 15. The spherical cavity had a radius of 1 cm. 


Table 15. Eigenvalues of Spherical Cavity 
With Radius of 1 cm 


Mode 

k Q , cm 1 



Reference 16 

TMoio 

2.744 

2.799 

2.799 

TMoio 

2.744 

2.802 

2.802 

TMqio 

2.744 

2.807 

2.811 

TM 021 

3.870 

3.961 

3.948 

TM121, even 

3.870 

3.976 

3.986 

TM121, odd 

3.870 

3.986 

3.994 

TM221, even 

3.870 

3.994 

4.038 

TM221, odd 

3.870 

3.998 

4.048 

TEO 011 

4.493 

4.421 

4.433 

TE111, even 

4.493 

4.478 

4.472 

TE111, odd 

4.493 

4.501 

4.549 


3 A. 6. Summary 

In section 3. 1 , vector finite elements introduced in 
section 2.2 are extended to solve three-dimensional 
eigenvalue problems by using tangential edge basis 
functions for tetrahedral elements. Spurious solu- 
tions are completely avoided because of the diver- 
gence free nature of the edge elements. Sparse matrix 
eigenvalue solvers are used to take advantage of the 
sparsity and symmetry of the finite element matrices; 
this results in considerable savings in computer mem- 
ory and computational time. The numerical results 
presented for cavities with different shapes prove the 
validity of the analysis and accuracy of the computer 
codes presented in this section. 

4. Concluding Remarks 

A thorough formulation of FEM for various eigen- 
value problems in electromagnetism has been pre- 
sented. The use of recently developed edge basis 
functions for vector finite elements has been demon- 
strated for two-dimensional and three-dimensional 
problems. Triangular elements for two-dimensional 
problems and tetrahedral elements for three- 
dimensional problems are used to model complex geo- 
metrical shapes because of their ability to represent 
such shapes accurately. 

Implementation of conventional nodal-based 
scalar formulation for two-dimensional homogeneous 
problems is demonstrated in section 2.1. The com- 
puter code developed can be used to calculate the 
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eigenvalues and field intensity patterns of any ar- 
bitrarily shaped waveguide filled with homogeneous 
material. The numerical results for various wave- 
guides and the field intensity patterns for various 
modes have been presented in section 2.1 which show 
the validity of the analysis and computer code. Sim- 
ple waveguide shapes have been chosen for demon- 
stration purposes because their field distributions are 
well known. 

In section 2.2, the two-dimensional edge basis 
functions for triangular elements are introduced to 
model transverse vector fields in waveguides. For 
waveguides or microwave circuit design problems in 
general, knowledge of a propagation constant at a 
given frequency is desirable. A step-by-step formu- 
lation is done to determine either the wave number 
or the propagation constant, if one of them is spec- 
ified. Numerical examples are presented to validate 
the analysis and the computer codes. 

Since real problems involve three-dimensional ge- 
ometries, a formulation for the calculation of eigen- 


values for such geometries is presented in section 3. 
The three-dimensional edge basis functions are in- 
troduced for tetrahedral elements. Numerical exam- 
ples for various geometries are presented. Compari- 
son of the numerical data with the available data in 
the literature shows the validity and accuracy of the 
analysis. 

For the three-dimensional problems, the number 
of variables increases drastically compared with those 
for the two-dimensional problems. In the three- 
dimensional computer codes, sparsity of FEM matri- 
ces have been exploited by storing only the nonzero 
elements and symmetry is utilized by storing only 
either the upper or lower part of the matrices. The 
sparse matrix eigenvalue solvers are used to efficiently 
solve the FEM equations. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
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Appendix 

Computer Programs 

Computer programs were written to implement 
the analysis presented in this report. All the pro- 
grams are written in FORTRAN language. These 
programs take *.MOD file containing the meshing 
information from COSMOS/M (ref. 4). Also these 
computer programs make use of the optimized library 
routines of EISPACK (refs. 7 and 8) and VECLIB 
(ref. 18) on the CONVEX computer. 

HELM 10 is a two-dimensional finite element pro- 
gram to calculate the eigenvalues of homogeneously 
filled waveguides. This program implements the 
analysis described in section 2.1. 

HELMVEC is a two-dimensional finite element 
program to calculate the eigenvalues of homo- 
geneously filled waveguides. Unlike HELM10, this 
program uses the vector basis functions and imple- 
ments the analysis presented in section 2.2.1. 

HELMVEC 1 is a two-dimensional finite element 
program to calculate the eigenvalues of inhomo- 
geneously filled waveguides. This program imple- 
ments the three-component vector field formulation 
discussed in section 2.2.2 and makes use of edge basis 
functions for transverse fields and scalar basis func- 
tions associated with nodes for longitudinal fields. 


HELMVEC2 is a two-dimensional finite element 
program for finding out the wave number of an in- 
homogeneously filled waveguide when the propaga- 
tion constant is specified. This program implements 
the formulation presented in section 2.2.3. 

HELMVEC3 is a two-dimensional finite element 
program for determination of propagation constant 
of an inhomogeneously filled waveguide at any fre- 
quency. This program implements the formulation 
given in section 2.2.4. 

FEM3D0 and FEM3D1 are three-dimensional fi- 
nite element programs to calculate the eigenvalues 
of inhomogeneously filled cavities. These programs 
implement the formulation presented in section 3.1. 
FEM3D0 uses the EISPACK (refs. 7 and 8) routines, 
whereas FEM3D1 exploits the sparsity and symme- 
try of FEM matrices and uses VECLIB (ref. 19) 
routines. 

These computer programs are available on request 
from 

Information and Electromagnetic Technology 
Division 

Electromagnetic Research Branch 

M.S. 490 

NASA Langley Research Center 

Hampton VA 23681-0001 


26 


References 

1. Silvester, P. P.; and Ferrari, R. L.: Finite Elements for 
Electrical Engineers. Second ed., Cambridge Univ. Press, 
1990. 

2. Rahman, B. M. Azizur; and Davies, J. Brian: Penalty 
Function Improvement of Waveguide Solution by Finite 
Elements. IEEE Trans. Microw. Theory & Techni ., 
vol. MTT-32, no. 8, Aug. 1984, pp. 922-928. 

3. Bossavit, Alain: Simplicial Finite Elements for Scattering 
Problems in Electromagnetism. Comput. Methods Appl. 
Mech . & Eng., vol. 76, 1989, pp. 299^316. 

4. COSMOS/M — User Guide. Volume 1 — Preprocessing, 
Analysis and Postprocessing Interface, Version 1.70. 
Structural Research and Analysis Corp., May 1993. 

5. I-DEAS Relational Data Base. I-DEAS V Prod. Bui., 
Structural Dynamics Research Corp., Jan. 1990. 

6. Silvester, P.: Construction of Triangular Finite Element 
Universal Matrices. Int. J. Numer . Methods Eng . vol. 12, 
no. 2, 1978, pp. 237-244. 

7. Smith, B. T.; Boyle, J. M.; Dongarra, J. J.; Garbow, 

B. S.; Ikebe, Y.; Klema, V. C.; and Moler, C. B.: Matrix 
Eigensystem Routines — EISPACK Guide. Second Ed., 
Springer- Verlag, 1976. 

8. Garbow, B. S.; Boyle, J. M.; Dongarra, J. J.; and Moler, 

C. B.: Matrix Eigensystem Routines — EISPACK Guide 
Extension. Springer- Verlag, 1977. 

9. Harrington, Roger F.; Time-Harmonic Electromagnetic 
Fields . McGraw-Hill Book Co., Inc., 1961. 

10. Marcuvitz, N., ed.: Waveguide Handbook. McGraw-Hill 
Book Co., Inc., 1951. 


11. Koshiba, Masanori; Hayata, Kazuya; and Suzuki, Michio: 
Improved Finite-Element Formulation in Terms of the 
Magnetic Field Vector for Dielectric Waveguides. IEEE 
Trans . Microw. Theory & Techni ., vol. MTT-33, no. 3, 
Mar. 1985, pp. 227-233. 

12. Whitney, Hassler: Geometric Integration Theory. 

Princeton Univ. Press, 1957. 

13. Lee, J.-F.; Sun, D. K.; and Cendes, Z. J.: Full-Wave 
Analysis of Dielectric Waveguides Using Tangential Vector 
Finite Elements. IEEE Trans. Microw. Theory & Techni ., 
vol. 39, no. 8, Aug. 1991, pp. 1262-1271. 

14. Reddy, J. N.: An Introduction to the Finite Element 
Method. McGraw-Hill Book Co., Inc., 1984. 

15. Hayata, Kazuya; Koshiba, Masanori; Eguchi, Masashi; 
and Suzuki, Michio: Vectorial Finite-Element Method 
Without Any Spurious Solutions for Dielectric Wave- 
guiding Problems Using Transverse Magnetic-Field Com- 
ponent. IEEE Trans. Microw. Theory & Techni ., 
vol. MTT-34, no. 11, Nov. 1986, pp. 1120-1124. 

16. Lee, Jin-Fa; and Mittra, Raj: A Note on the Applica- 
tion of Edge-Elements for Modeling Three-Dimensional 
Inhomogeneously-Filled Cavities. IEEE Trans. Microw. 
Theory & Techni vol. 40, no. 9, Sept. 1992, 
pp. 1767-1773. 

17. Chatterjee, A.; Jin, J. M.; and Volakis, J. L.: Computa- 
tion of Cavity Resonances Using Edge-Based Finite Ele- 
ments. IEEE Trans. Microw. Theory & Techni ., vol. 40, 
no. 11, Nov. 1992, pp. 2106-2108. 

18. Zienkiewicz, O. C,: The Finite Element Method in Engi- 
neering Science. McGraw-Hill Book Co., Inc., 1971. 

19. CONVEX VECLIB—User’s Guide , Seventh ed., 
CONVEX Computer Corp., 1993. 


27 



TTT F 


REPORT DOCUMENTATION PAGE 

Form Approved 
OMB No. 0704-0188 

Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwprk Reduction Project (0704-0188), Washington, DC 20503 

1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 

December 1994 Technical Paper 

4. TITLE AND SUBTITLE 

Finite Element Method for Eigenvalue Problems in Electromagnetics 

5. FUNDING NUMBERS 

WU 505-64-52-04 

6. AUTHOR(S) 

C. J. Reddy, Manohar D. Deshpande, C. R. Cockrell, and Fred B. Beck 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

NASA Langley Research Center 
Hampton, VA 23681-0001 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

L- 17392 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 

10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

NASA TP-3485 

11. SUPPLEMENTARY NOTES 

Reddy: NRC-NASA Resident Research Associate, Langley Research Center, Hampton, VA; Deshpande: 
ViGYAN, Inc., Hampton, VA; Cockrell and Beck: Langley Research Center, Hampton, VA. 

12a. DISTRIBUTION/AVAILABILITY STATEMENT 1 12b. DISTRIBUTION CODE 


Unclassified-Unlimited 
Subject Category 32 

Availability: NASA CAST (301) 621-0390 


13. ABSTRACT (Maximum 200 words) 

Finite element method (FEM) has been a very powerful tool to solve many complex problems in electromag- 
netics. The goal of the current research at the Langley Research Center is to develop a combined FEM/method 
of moments approach to three-dimensional scattering/radiation problem for objects with arbitrary shape and 
filled with complex materials. As a first step toward that goal, an exercise is taken to establish the power of 
FEM, through closed boundary problems. This paper demonstrates the development of FEM tools for two- 
and three-dimensional eigenvalue problems in electromagnetics. In section 2, both the scalar and vector finite 
elements have been used for various waveguide problems to demonstrate the flexibility of FEM. In section 3, 
vector finite element method has been extended to three-dimensional eigenvalue problems. 


14. SUBJECT TERMS 

Finite element; Waveguide; Eigenvalues; Computational method; Electromagnetic 


15. 


NUMBER OF PAGES 

36 


16. 


PRICE CODE 

AQ3 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 


NSN 7540-01-280-5500 


18. SECURITY CLASSIFICATION! 
OF THIS PAGE 

Unclassified 


19. SECURITY CLASSIFICATION! 
OF ABSTRACT 

Unclassified 


20. LIMITATION 
OF ABSTRACT 


Standard Form 298(Rev. 2-89) 

Prescribed by ANSI Std Z39-18 
298-102 





















FT i 



